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APPENDIX A 
% BEGINNING OF PSEUDO CODE 

5 % compute scale factor A, and time constants a, b from physical system 
% parameters 

A = Vmax * Kt / (Re * Rm + Kt * Kb) * 1 * k; 

10 pi = 1/Jm/Ie * (-Ie * Rm - Re * Jm + sqrt(Ie"2 * Rm"2 - 2 * Re * Rm * Ie 

* Jm + Re"2 * Jm"2 - 4 * Kt * Kb * Ie * Jm)) / 2; 

p2 = 1/Jm/Ie * (-Ie * Rm - Re * Jm - sqrt(Ie A 2 * Rm A 2 - 2 * Re * Rm * Ie 

* Jm + Re^2 * Jm*2 - 4 * Kt * Kb * Ie * Jm)) / 2; 

15 a = max(-pl,-p2) 

b = min(-pl,-p2) 

% make initial guesses for step durations 

20 etl = 1; 

et2 = .005; 
et3 = 1; 

% set maximum iteration count 

25 

Nmax = 1000; 

forj = l:Nmax 

% save old values of step time intervals 
30 et3old = et3; 
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et2old = et2; 
etlold = etl; 

% iterate for switch times using fixed voltage level Vmax 

et3 = -log(1.0 / 2.0 - exp(-etl * a) / 2 + exp(-et2 * a)) / a; 

et2 = 1/b * log(2.0) + 3 * et3 - 1/b * log(2 * exp(l/A * b * X) * 

exp(et3 * b) - sqrt( 4.0) * sqrt(exp(l/A * b * X)) * exp(et3 * b) 
* sqrt(exp(l/A * b * X)+exp(et3 * b)"2 - 2 * exp(et3 * b))); 

etl = - (-2 * A * et2 + 2 * A * et3 - X) / A; 

if norm([et3old - et3 et2old - et2 etlold - etl], inf) < = eps * 2 
break 

end 

if j = = Nmax 

error(['error - failure to converge after ', num2str(Nmax),' 
iterations']) 

end 

end 

% round up pulse duration to nearest sample interval, 

% convert to intervals between steps to make sure that voltage 

% requirements will not increase (beyond Vmax). 

dtl=ceil((etl - et2) / dt) * dt; 
dt2=ceil((et2 - et3) / dt) * dt; 
dt3=ceil((et3) / dt) * dt; 

etl23 = [etl, et2, et3] 

% convert back to total step duration. 



etl = dtl + dt2 + dt3; 
et2 = dt2 + dt3; 
et3 = dt3; 

5 % In the following, the original constraints equations involving XF1, XF2, 

% and XF3 have been modified to include a variable voltage level applied at 
% each step (instead of the fixed maximum (+/-) Vmax). 

% The original equations for XF1, XF2, and XF3 follow: 

10 % XF.CU) = VoF.Ct^ - to) - 2V <F it end - 1 j) + 2V <F j[t end - 1 } 

% XF 2 (t end ) = VoF^ - g - 2V 0 F 2 (t end - 1,) + 2V 0 F 1 (t end - 1 2 ) 

% XF 3 (t end ) = V 0 F 3 (t tend - g - 2V 0 F 2 (t end - 1,) + 2V 0 F,(t end - 1 2 ) 

% And the modified equation including adjustable relative levels of voltage 
15 % LI, L2 and L3 are: 

% XF,(t end ) = L.VoFAend - g - L2V 0 F,(t end - g + L 3 V<F£ end - 1) 
% XF 2 (t end ) = L.VoF^ - g - L 2 Vf 2(t end - ti) + L 3 V^{t end - t^ 
% XF 3 (t end ) = L 1 V 0 F 3 (t tend - g - LjVoF^ - 1,) + I^VoF,^ - 1 2 ) 

20 % And the corresponding constraint equations are: 
% XF,(t end ) = Finalpos 
% XF 2 (t end ) = 0 
% XF 3 (t end ) = 0 

25 % Where all of the times indicated have discrete values, e.g. corresponding to 

% the controller update rate. 

% It should be noted that after the digital switch times are fixed, the constraint 
% equations derived from the equations above form a linear set of equations in 
30 % the unknown relative voltage levels LI, L2 and L3 and any standard linear 
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% method can be used to solve for the relative voltage levels. In the equations 
% for (LI, L2 and L3) that follow, the solution was obtained by algebraic 
% means (and are not particularly compact.) 



5 % compute new relative voltage step levels 

% LI, L2 and L3 are nominally assigned to "1", "-2" and "+2", respectively 



si = X * (exp(-et3 * b) * exp(-et2 * a) + exp(-et3 * a) + exp(-et2 *b ) - 
exp(-et2 *b) * exp(-et3 * a) - exp(-et2 * a) - exp(-et3 *b )); 

10 s2 = 1 / (et2 * exp(-etl * b) * exp(-et3 * a) + exp(-et2 * b) * et3 * 

exp(-etl * a)- et2 * exp(-et3 * a) - et2 * exp(-etl * b) -et3 * 
exp(-etl * a) - exp(-et2 * b) * et3+exp(-et3 * b) * etl * exp(-et2 * a) 
+ exp(-et3 * a) * etl + exp(-et2 * b) * etl - exp(-et2 * b) * etl * 
exp(-et3 * a) - et3 * exp(-etl * b) * exp(-et2 * a) - exp(-et2 * a) * etl - 

15 exp(-et3 * b) * etl - exp(-et3 * b) * et2 * exp(-etl * a) + et3 * 

exp(-etl * b) + et2 * exp(-etl * a) + exp(-et3 * b) * et2 + et3 * 
exp(-et2 * a)) / A; 



LI = si * s2; 

20 

si = 1 / (et2 * exp(-etl * b) * exp(-et3 * a) + exp(-et2 * b) * et3 * 
exp(-etl * a) - et2 * exp(-et3 * a) - et2 * exp(-etl * b) - et3 * 
exp(-etl * a) - exp(-et2 * b) * et3 + exp(-et3 * b) * etl * 
exp(-et2 * a) + exp(-et3 * a) * etl+exp(-et2 * b) * etl - 

25 exp(-et2 * b) * etl * exp(-et3 * a) - et3 * exp(-etl * b) * 

exp(-et2 * a) - exp(-et2 * a) * etl - exp(-et3 * b) * etl - exp(-et3 * b) * 
et2 * exp(-etl * a) + et3 * exp(-etl * b) + et2 * exp(-etl * a) + 
exp(-et3 *b ) * et2 + et3 * exp(-et2 * a)) * X; 
s2 = (exp(-et2 * b) * exp(-etl * a) - exp(-etl * a) - exp(-et2 * b) - 

30 exp(-etl * b) * exp(-et2 * a) + exp(-etl * b) + exp(-et2 * a)) / A; 
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L3 = sl*s2; 

si = exp(-etl * a) - exp(-et3 * a) + exp(-et3 * b) - exp(-etl * b) - 
exp(-et3 * b) * exp(-etl * a) + exp(-etl * b) * exp(-et3 * a); 
5 s2 = X / (et2 * exp(-etl * b) * exp(-et3 * a) + exp(-et2 * b) * et3 * 
exp(-etl * a) - et2 * exp(-et3 * a) - et2 * exp(-etl * b) - et3 * 
exp(-etl * a) - exp(-et2 * b) * et3 + exp(-et3 * b) * etl * exp(-et2 * a) 
+ exp(-et3 * a) * etl + exp(-et2 * b) * etl - exp(-et2 * b) * etl * 
exp(-et3 * a) - et3 * exp(-etl * b) * exp(-et2 *a ) - exp(-et2 * a) * etl- 
10 exp(-et3 * b) * etl - exp(-et3 * b) * et2 * exp(-etl * a) + et3 * 

exp(-etl * b) + et2 * exp(-etl * a) + exp(-et3 * b) * et2 + et3 * 
exp(-et2 * a)) / A; 

L2 = si * s2; 

15 

% convert accumulated voltage steps to sequential voltage level 

VI = Vmax * (LI); 

V2 = Vmax * (LI + L2); 

V3 = Vmax * (LI + L2 + L3); 

20 

% END OF PSEUDO CODE 
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